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Abstract 

The contributions a nj to the O(a^) massive operator matrix elements describing the 
heavy flavor Wilson coefficients in the limit Q 2 3> m 2 are computed for the structure 
function F2(x,Q 2 ) and transversity for general values of the Mellin variable N. Here, 
for two matrix elements, A^q(N) and A qg: Q(N), the complete result is obtained. A first 
independent computation of the contributions to the 3-loop anomalous dimensions j qg (N), 
7qq(N), and "fqq 5 '^ TR \N) is given. In the computation advanced summation technologies 
for nested sums over products of hyper geometric terms with harmonic sums have been 
used. For intermediary results generalized harmonic sums occur, while the final results can 
be expressed by nested harmonic sums only. 



1 Introduction 



The heavy flavor corrections to deep-inelastic structure functions amount to large contributions 
at lower values of the Bjorken variable x. Currently they are known in semi-analytic form to 
2-loop (NLO) order [TJ. The present accuracy of the deep-inelastic data reaches the order of 
1% [2], which requires the next-to-next-to-leading order (NNLO) corrections for precision deter- 
minations of both the strong coupling constant a s (M§) and the parton distribution functions [3], 
as well as the detailed understanding of the heavy flavor production cross sections in lepton- 
nucleon scattering [4]. The precise knowledge of these quantities is of central importance for 
the interpretation of the physics results at the Large Hadron Collider, LHC, [5]. In the region 
Q 2 3> m 2 , with Q 2 = —q 2 , with q the space-like 4-momentum transfer and m the heavy quark 
mass, the power corrections 0((m 2 /Q 2 ) k ), k > 1 to the heavy quark structure functions become 
very small. For the structure function F 2 (x,Q 2 ) the logarithmic and constant contributions are 
sufficient at the 1%- level to describe the complete result for Q 2 /m 2 > 10, a region which does 
well compare to the deep-inelastic region at HERA in which the twist- 2 contributions dominate, 
cf. 0.0 In this limit the Wilson coefficients with rif massless and one massive quark factorize 
into massive operator matrix elements (OMEs) and the massless Wilson coefficients, as has been 
shown in Ref. [8]. The former quantities are process independent, while the latter depend on 
the respective scattering process. The massless Wilson coefficients for the structure function 
F 2 (x, Q 2 ) are known to 3-loop order, [9]. 

For fixed Mellin moments N a series of moments up to iV = 10. ..14, depending on the 
respective transition, have been calculated for all the OMEs at 3-loop order contributing to the 
structure function F 2 (oc,Q 2 ) and those needed to establish a variable flavor scheme description 
at 0(otg) in Ref. [10]. □ There also the complete renormalization of the matrix elements has 
been derived. In this computation the massive OMEs for given total spin N were mapped onto 
massive tadpoles which were computed using MAT AD, [12]. For general values of N the 2-loop 
OMEs, up to 0(e), have been calculated in Refs. [K5HT5] . All the logarithmic contributions 
to the massive OMEs are known [16|[T7] for general values of N. They depend on the 3-loop 
anomalous dimensions [IEJ[19]. F° r the structure function Fl(x, Q 2 ) the asymptotic heavy flavor 
Wilson coefficients at 0(a%) were calculated in [20]. They become, however, effective only at 
much higher scales of Q 2 compared to the case of F 2 (x, Q 2 ). 

In the present paper the 0(a^) contributions oc ripTpCp^A are computed for all massive 
operator matrix elements contributing to the structure function F 2 (x, Q 2 ) at general values of the 
Mellin variable N in the fixed flavor number scheme, as well as the corresponding contributions 
to transversity. This scheme has to be considered as the genuine scheme in quantum field 
theoretic calculations since the initial states, the twist-2 massless partons can, at least to a good 
approximation, be considered as LSZ-states. This is not the case for heavy quark states, which 
have a finite lifetime. For two OMEs, A™ Q (N) and A qg>Q (N), the complete result is obtained. 
In the present computation the Feynman parameter integrals are computed directly. They can 
be represented in terms of generalized hypergeometric functions [22J and sums thereof prior the 
expansion in the dimensional variable e = D — 4, cf. [23j[2l]. Finally, they are represented in 
terms of nested sums over products of hypergeometric terms and harmonic sums, which can be 
calculated using modern summation techniques [251 ES] that are based on a refined difference 
field of [27] and that generalize the summation paradigms presented in [28] to mult i- summation. 

x For higher order corrections to the gluonic contributions in the threshold region, cf. [7]. 
2 For the corresponding contributions in case of transversity see 

3 From the representations obtained in the fixed flavor number scheme, variable flavor number schemes may 
be defined under specific conditions [T0J[T3] observing the correct matching scales |21) . 
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During this computation the results can be expressed in terms of nested harmonic sums [29,30j. 
In intermediary steps of the calculation generalized harmonic sums, [3TJ|32] , cf . also |33j , appear 
which finally cancel. 

The paper is organized as follows. In Section [2] we summarize the basic formalism. The 
results for the constant part of the O(a^) n/-contributions to the massive OMEs Aq 9 (N), Aq(N) 
Aqq, Q (N), A qgtQ (N), A p q l Q {N), and A^ R {N), cf. HDUm, are presented in Section! The single 
pole terms in e allow to derive the terms oc rif of the 3-loop anomalous dimensions for general 
values of N. They are compared to the results in Refs. [18l[T9],[3l] and are obtained in a first 
independent calculation for j gg (N), 7^ s (iV), 7^ 9 S,TR (A^), in Section HI Section [5] contains the 
conclusions. Some technical details of the calculation are given in the Appendix. 



2 The heavy flavor Wilson coefficients in the asymptotic 
region 

The heavy flavor contributions to the structure functions F^,l)(x, Q 2 ) with rif massless and one 
heavy flavor are given by, pH] : 



F® Q L) (x,n f + l,Q 2 ,m 2 ) = ^ e£ jlrj^j (x, n f + 1, , ^ ) ® f k (x, n\ n f ) + /~(z, ^, n f ) 



l r 
+— 



L 



PS 

9,(2,L) 



Q 2 m 2 

,n,f + 1, —j, —j ) ® T,(x,fi 2 ,nf) + L 

12 ^,2 



+4 



3>(2,i) 



g2 m 2 



PS f Q ™ A / 2 \ c ( 



Q 2 ^ 
9 5 9~ 



G(x,n 2 ,n f ) 



I) 



The different Wilson coefficients are denoted by Li, Hi in case the photon couples to a light (L) 
or the heavy (H) quark, for the flavor non-singlet (NS), pure-singlet (PS), and singlet (S) cases. 
Here, <g> is the Mellin convolution, 



[A <g> B](x) 



1 rl 



JO 



dx\dx2 S(x — X1X2) A(xi) B (X2) 



(2) 



with boundaries for the Wilson coefficients [x(l + Am 2 /Q 2 ), 1], e k the light and cq the heavy 
quark charges, /i 2 denotes the factorization scale, and £ and G are the quark, antiquark, 

flavor singlet and gluon momentum distribution functions, with 



E(x, /i 2 , n f ) = [fk(x, /i 2 , n f ) + fe(x, /z 2 , n/)] . 



(3) 



k=l 



For Q 2 ^> m 2 the massive Wilson coefficients can be expressed in terms of the renormalized 
massive OMEs A^ and the massless Wilson coefficients Cj. To 0(0%) they read (a s = a a /(47r)), 

cf. nm : 



L NS 



g,(2,i) 



+ a* 



PS 



Uf + l) 



A^f + ^ + C^(nf) 

A^ 5 (n f + % + ^f(n, + 1)^>/ + 1) + C^(n f ) 
42q>/ + 1 ) 5 * + A ?1qM n Al, L) (n f + 1) + n f df}( S L) (n f ) 
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L l(2,L)( n f + 1) 



(i) 



V(2) 



(2) 



(3) 



fn f + 1) <5 2 



n f + L)+n f C gM {nf) 



+ 4&(n, + 1) n f C 



H q S (2M"J 



n f + l) 



l(2),PS, 
V(3),PS , 



.4 



(3), PS 



(71/ + 1) <5 2 



(i) 



+^ PS (n/ + l)C}^ ) (n / + l) 



(1),NS. 



S,(2,i) 



(n/ + 1) 



(1) 

9,(2,L)' 



AgVn/ + 1) 5 2 



+ + 1) ^(n, + 1) + A gg [ Q (n f + 1) Cj^n, + 1) 



(i) 



+ CS,L)(% + 1)J + < [ 4l(n f + 1)6 2 + A$ g (n f + 1) C™*(n, + 1) 
+ Afl Q {n, + 1) c£l L) (n f + 1) + Ag](n/ + l){c$S(n, + 1) 
+ + 1)} + A« Q (n/ + 1) 0^,(71/ + 1) + <7$,, L) (n, + 1) 



(4) 

with 6 2 = 0(1) for F L (F 2 ) and f(n f ) = f(n f + 1) - /(n/), f(n f ) = f(n f )/n f . 

The renormalized massive OMEs depend on the ratio m 2 //i 2 , while the scale ratio in the 
massless Wilson coefficients is fi 2 /Q 2 . The latter are pure functions of the momentum fraction 
z, or the Mellin variable N, if one sets fi 2 = Q 2 . The mass dependence of the heavy flavor Wilson 
coefficients in the asymptotic region derives from the unrenormalized massive OMEs 



(3),3 



3% 



i(3),2 



(3),1 



i(3),0 

'ij 



applying mass, coupling constant, and operator-renormalization, as well as mass factorization, 
cf. Ref. [Eg. 

The renormalized massive OMEs obey then the general structure 

,2\ / m 2\ /™2\ /™2\ 

(6) 



A?) ' '"' 



Q 2 



ar ,3g) +0 -^g) +a g^(^) + ^ 



The subsequent calculations will be performed in the MS scheme. For other scheme choices see 
Ref. [10J. Therefore the strong coupling constant is obtained as the perturbative solution of the 
equation 



da s {fi 2 
rfln(yU 2 



(7) 



to 3-loop order, where /3 k are the expansion coefficients of the QCD /3-function and /i 2 denotes 
the renormalization scale. For simplicity we identify the factorization and renormalization scales 
in the following. 



3 The Massive Operator Matrix Elements 

The operator matrix elements oc n/ for both F 2 (x,Q 2 ) and transversity are obtained by the 
massive two-loop graphs inserting a further massless fermion line and new planar three-loop 
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topologies, cf. [SSIEE], as well as 3-loop graphs containing bubble topologies with operator inser- 
tions linked either linked to massive or massless fermion lines, cf . [10] . The calculation was carried 
out in Feynman-gauge using FORM [37] and MAPLE-codes, and applied the package color [38] 
for the color algebra. As in earlier cases [15] we computed the Feynman parameter-integrals di- 
rectly, without applying the integrating-by-parts method [39]. The corresponding integrals can 
be mapped onto sums over generalized hypergeometric functions prior the e-expansion, which 
allow to obtain the Laurent series in e. Finally, up to three-fold nested sums over hypergeo- 
metric expressions, equipped with harmonic sums, have to be performed, for which the package 
Sigma [25], constructing difference and product fields, was applied and extended. Some details 
of the computation are presented in Appendix A. 

The massive OMEs (N) are finally obtained as rational functions of the Mellin variable 
N, multiple zeta values |40j, and nested harmonic sums [291130] . The latter are defined recursively 
by 

S b ,a(N) = (sig ^ 6))fc ^(fc), S 9 (N) = 1 . (8) 
fc=i 

As a short-hand notation we use Sg(N) = S3. In representing the results, the algebraic relations 
of the nested harmonic sums |4T] are applied. In the following we present the constant con- 
tributions to the unrenormalized OMEs (jSJ) as genuine quantities, to allow for different scheme 
choices, cf. Ref. HU1. 



3.1 Operator Matrix Elements contributing to F 2 (x,Q 2 ) 

The 0(nf) contribution to the unrenormalized OME Aq (s,N), clq' , reads : 

fig' = n f 7fC A i MNW + tyN ' ^ \ mS ~W ~ 78 ^,i,i ~ 905-3,! + 72S^ 2 - 6^1 



+ 2 



2 



-1085_ 2 ,iSi + 425 2 ,i5i - 6S_ 4 + 905_ 3 ^i + HSS 3 S 1 + 120S 4 + ISS. 2 S 2 + MS- 2 S( 
+33S 2 Sl + 15S 2 2 + 2Sf + 185_ 2 C 2 + 9S 2 ( 2 + 95^2 - 42Si( 3 
5iV 4 + 14iV 3 + 53iV 2 + 82iV + 20 



+32 



27iV(iV + l) 2 (iV + 2) 2 
64(5iV 4 + lliV 3 + 50iV 2 + 85iV + 20) 

27iV(iV + l) 2 (iV + 2) 2 
16(40iV 4 + 151iV 3 + 544iV 2 + 779N + 214) 

27N(N +1) 2 (N + 2) 2 
32(65iV 6 + 429iV 5 + 1155iV 4 + 725iV 3 + 370iV 2 + 496iV + 648 



QS— 2 1 — 55—3 — QS— 2 S\ 
S 2 ,i 

n 

S 2 Si 

I- 4Pfi N 4- f\dR) 

s* 



2 



81(N - l)iV 2 (iV + l) 2 (iV + 2) 
16(20iV 4 + 107iV 3 + 344iV 2 + 439iV + 134) 3 Q x (N) 



81iV(iV + l) 2 (iV + 2) 2 1 81(JV- 1)N 3 (N + l) 3 (iV + 2) 3 

32(47iV 6 + 278iV 5 + 1257iV 4 + 2552iV 3 + 1794iV 2 + 284iV + 448) 
+ 81iV(iV + l) 3 (iV + 2) 3 ^ 



4 In Ref. [10] we have kept the gauge parameter for part of the moments and found gauge independence. In 
the present calculation we have compared the results also on the basis of diagrams with the moments obtained 
there. 
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8(22iV 6 + 271iV 5 + 2355iV 4 + 6430iV 3 + 6816iV 2 + 3172iV + 1256) 2 

<->i 



3 



81iV(iV + l) 3 (iV + 2) 
Q 2 (N) U8(N 2 + N + l)(N 2 + N + 2) 

L^l + nf at nAraar , i\our , r>\2 ^ 3 



243(iV- l)iV 2 (iV+ l) 4 (iV + 2) 4 1 9(iV- l)A^ 2 (iV+ l) 2 (iV + 2) 
16(5iV 4 + 20iV 3 + 59iV 2 + 76iV + 20) Q 3 (N) 



9iV(iV + l) 2 (iV + 2) 2 ^ z 9(iV- l)iV 3 (iV + l) 3 (iV + 2) 3 

Q 4 (AQ 

243(iV - l)iV 5 (iV + l) 5 (iV + 2) 5 

+n/r^ 1 27 ^7iXJV + 2) I 14452 ' 1 - 1 - 72 ^ 3 ' 1 " 7252 ^ + 4854 - 165351 

-24S 2 - 12^ 2 - 2Sf - 9^ 2 C 2 + 42S&] + 32^ + ^ + ^ \ ^ hs^ + Sf] 
2 i i iS2 isjj 81iV 2 (iV + l)(7V + 2) L : J 

128(iV 2 - 3N - 2) Q 5 (iV) 



3iV 2 (iV + l)(iV + 2) 2,1 ~~ 81(iV - 1)N 3 {N + 1) 3 {N + 2) 2 3 

Q 6 (N) 32(10iV 4 + 185iV 3 + 789iV 2 + 521iV + 141) 2 

+ 27(iV - l)iV 4 (iV + l) 4 (iV + 2) 3 2 81iV 2 (iV + l) 2 (iV + 2) 1 

16(230iV 5 - 924iV 4 - 5165iV 3 - 7454iV 2 - 10217iV - 2670) 

243iV 2 (iV + l)3(iV + 2) 1 
16(5iV 3 + lliV 2 + 28iV + 12) Q 7 (N) 



9iV 2 (iV + l)(iV + 2) b 9(7V- l)iV 3 (iV + l) 3 (iV + 2) 

Qs(iV) g 9 (iV) 
rs2 + 



9(iV- l)iV 4 (iV+ l) 4 (iV + 2) 3 ^ 243(iV- l)iV 6 (iV + l) 6 (iV + 2) E 
with the polynomials 

Qi{N) = 32N 9 - 936N 8 + 6448iV 7 + 55208iV 6 + 126160iV 5 + 61760iV 4 

-53152iV 3 - 250247V 2 - 32256iV - 13824 , 
Q 2 (N) = +7856 TV 10 + 84672 iV 9 + 377648/V 8 + 985568 iV 7 + 1395456/V 6 

+470688/V 5 - 1183712iV 4 - 1180224/V 3 - 182528/V 2 - 42752/V 

+ 13824 , 

Q 3 (N) = QON 9 + 360iV 8 + 584iV 7 - 128iV 6 - 2004iV 5 - 2440iV 4 - 976/V 3 

-192iV 2 + 896iV + 384 , 
Q 4 (N) = 28776iV 15 + 356112Ar 14 + 1896088Ar 13 + 5538320A^ 12 + 9112264A^ 11 

+6793968iV 10 - 3019528iV 9 - 11879520iV 8 - 11673088/V 7 

-6450992iV 6 - 3726976iV 5 - 2248128/V 4 - 183296/V 3 + 268032iV 2 

+147456iV + 27648 , 
Q 5 (JV) = +464A^ 8 -15616iV 7 -38112A^ 6 + 27776iV 5 + 146064A^ 4 + 119552A^ 3 

+ 109312/V 2 + 86016iV + 62208 , 
Qq(N) = 456iV n + 4376iV 10 + 11328iV 9 - 3184/V 8 - 54552iV 7 - 111720iV 6 

-155376iV 5 - 251072iV 4 - 312192iV 3 - 222464iV 2 - 135936iV 

-41472 , 

Q 7 (N) = 1Q8N S + 672iV 7 + 784iV 6 - 3192iV 4 - 5600iV 3 - 7168iV 2 - 4480iV 



Qs(N) 
Q 9 {N) 



-2688 , 

90iy n + 630iV 10 + 1592iV 9 + 126(W 8 - 1934iV 7 - 8218iV 6 
-15524iV 5 - 23944N 4 - 26752iV 3 - 18400iV 2 - 11328iV - 3456 , 
15777iV 17 + 186525iV 16 + 879391iV 15 + 1874085AT 14 + 575913iV 13 
-5568833N 12 - 10465411iV 11 - 2970289iV 10 + 11884298N 9 
+ 12640320N 8 - 10343664iV 7 - 40750480N 6 - 55711424N 5 
-53947712iV 4 - 42534912N 3 - 23256576iV 2 - 7865856iV 
-1244160 , 



(15) 
(16) 
(17) 



and 



1 



i=i 



(18) 
(19) 



denotes the Riemann (^-function. 

The corresponding contribution to the pure singlet OME A^'^\e,N) is given by 



a 



PS,(3),0 
Qq 



n f T 2 F C F 



N 2 {1 + N) 2 (2 + N)(N - 1) 



(N 2 + N + 2f 



1760 208 16 3 

<->3 7T~ <->2<->l — — D 1 



27 



9 



27 



™q<> , 224 r \ | Qio W 

3 g S3 I f ^(1 + ^(2 + ^) 



32 



Qn(iV) 



A + 



32 



208 16 2 16 

Qi 2 (iV) 



with 



81 N 2 (l + N) 2 (2 + N) 2 243 N 3 (I + N) 3 (2 + NY 

Q 10 (N) = 8N 7 + 37iV 6 + Q8N 5 -lliV 4 - 8Q N 3 -56iV 2 -1047V -48 
Qu(N) = 25 N 10 + 176 N 9 + A17 N 8 + 30 N 7 - 20 N 6 + 18A8 N 5 

+2244 N 4 + 1648 N 3 + 3040iV 2 + 2112iV + 576 , 
Q 12 (N) = 158 iV 13 + 1663 A^ 12 + 7714 iV 11 + 23003 A^ 10 + 56186 iV 9 
+89880 N s + 59452 N 7 - 8896 N 6 - 12856 N 5 - 24944 N 4 
-84608 N 3 - 77952 N 2 - 35712 N - 6912 , 

PS,(3), 



(20) 

(21) 
(22) 

(23) 



The second pure-singlet operator matrix element is A qq q (e, N). Its constant term reads : 



~PS,(3),0 
%q,Q 



n f T2 F C F ( 2 2 / 256 128 

N 2 (N-1)(2 + N)(1 + N) 2 \ {N +N + 2) l^T 3+ 9 2 ^ 



128 n , 32 224 



Qi 3 (N) 



64 



Qu{N) 



N(2 + N)(l + N) 
32 Q 15 (N) 



64 n 64 ^ 16 A 
— S 2 + — 5? + — C2 
27 27 1 9 s 



81 A^ 2 (2 + A^) 2 (l + A^) 2 243 N 3 (2 + N) 3 (l + N)'* 



(24) 



with 



Q 13 (N) = 16N 7 + 7AN 6 + 181 + 266 N 4 + 269 N 3 + 230 N 2 
+44 N - 24 , 



(25) 
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Qu(N) = 181 N 10 + 1352 N 9 + 4737 N 8 + 10101 N 7 

+14923 N 6 + 17085 N 5 + 14133 iV 4 + 5944 N 3 + 568 N 2 - 48 iV + 144 , (26) 
Q 15 {N) = 2074 N 13 + 21728 N 12 + 105173 iV 11 + 311482 N 10 + 636490 

+966828 + 1126568 N 7 + 968818 N 6 + 550813 N 5 

+ 169250 N 4 + 12104 N 3 - 3408 N 2 - 1008 N - 864 . (27) 
The constant term of the unrenormalized flavor non-singlet operator matrix element 

.ns,(3),o n T 2 C /64 448 32 320 

160 112 3N 2 + 3N + 2 640 
C2 ^1 ~ ; 7» ta at — C3 + t^T ^2 



27 ^ 27 (1 + N)N ™ 27 

4 3AT 4 + 6Ar 3 + 47Af 2 + 20Af- 12 55552 
+ 27 (1 + N) 2 N 2 C2_ ^29~ 1 

l 2 Qie(N) ] 
729 (1 + N^N* J ' 1 J 



where 



Q 16 (JV) = 1 1 75 1N 8 + 47004 N 7 + 93754 N 6 + 104364 N 5 + 55287 N 4 

+6256 N 3 - 2448 N 2 - 144 iV - 432 . (29) 

Finally, the constant contribution to A^q(s, N) at 0(rif) is : 



a (3),o _ n f T 



2 



N(N + l)(JV + 2) 

|cf [(iV 2 + ^ + 2) (1^ + j^S? + ls 2 S* + ^Ca^i + + \S 2 2 

,0, o, 40 o 56Q 17 (AQ x 

+8C 2 s 2 + y s 4 - 9(iV _ 1)iV2(iV + 1)2(iV + 2) C3j 

16(HW 3 + 13N 2 + 29N + 6)Sf 8(215N 4 + 481N 3 + 930N 2 + 748N + 120) 2 

8L/V + 81N(N + 1) 1 

16(HW 3 + 13N 2 + 29N + 6) / \ 32(40A^ 3 + 61N 2 + 89N + 6) 

27iV l 3 ^ 1 + 5251 J 81A7 S '- 

8(221A^ 4 + 515A^ 3 + 814AT 2 + 548A^ + 40) 4Q 18 (A0 

+ 27N(N + 1) 2 + 9(N - 1)N 3 {N + 1) 3 (A^ + 2) 2 ^ 2 

16Qi 9 (A0 8Q 20 (N) 



3 



243N(Af + l) 2 243(N- 1)N 5 (N + l) 5 (N + 2) 4 . 
+C A [(N 2 + AT + 2) (-±S* - ? C 2^ + jj^S? - fs 4 - iH?5 s>1 + |s 2>1>1 

H — zrS'iSi — — S01S1 — ~Sn — 5*_4 — — Cs5i — — C2S— 2 — — C2S0 

27 9 9 9 9 3 3 

448(N 2 + AT + 1) A \ + 32(5N 4 + 20N 3 + 41N 2 + A9N + 20) ^ 



:Ca) 



9(JV- l)iV(iV + l)(JV + 2) 81(iV + l)(iV + 2) 
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with 



x - SS 2 S 1 + US 2 ,i + sfj + 



128 
l27~ 



(5iV 2 + 8N + 10)5- 
8Q 22 (N) 



64 (5iV 4 + 38iV 3 + 59iV 2 + 31iV + 20) 
81(iV + l)(iV + 2) 

*Q2i(N) 



\1(N + l) 2 (iV + 2) : 



,D 1 



9(N -1)N 2 (N +1) 2 (N + 2) 2 
8Q 23 (N) 



32(i21iV 3 + 293iV 2 + AUN + 224) 

C2 777m~, ~ ^-2 



81(iV+ 1) 



\1{N + l) 2 (iV + 2) 2 

i6g 25 (iv) 



So 



ieg 24 (iv) 



243(iV - 1)N(N + 1) 3 (N + 2) 3 



5i 



243(iV - 1)N 4 (N + l) 4 (iV + 2) 1 



Qis(N) 

g 20 (iv) 



Q 22 (iV) 

Q 23 (iV) 
Q 24 (iV) 

Q 25 (iV) 



3iV 6 + 9N 5 — N 4 — 17N 3 - 38N 2 - 28N - 24 , 

18JV 11 + 126iV 10 + 365iV 9 + 630iV 8 + 652iV 7 + 626iV 6 + 1309iV 5 

+3170iV 4 + 4736iV 3 + 3584iV 2 + 2352iV + 864 , 

2507iV 5 + 8076iV 4 + 16120iV 3 + 18997iV 2 + 9898iV + 1344 , 

2322iV 17 + 30186iV 16 + 177047iV 15 + 627060iV 14 + 1509207iV 13 

+2623160iV 12 + 3436402A^ n + 3728602iV 10 + 4151281iV 9 

+5013306iV 8 + 5011065iV 7 + 3770902iV 6 + 3291500iV 5 + 3951272iV 4 

+3797616iV 3 + 2319264iV 2 + 862272iV + 155520 , 

2067V 6 + 1361iV 5 + 4134iV 4 + 7577iV 3 + 8394iV 2 + 4868iV + 1144 , 

6iV 9 + 36iV 8 + lliV 7 - 257iV 6 - 825iV 5 - 1375iV 4 - 1396iV 3 - 984iV 2 

-352iV - 48 , 

332iV 6 + 2537iV 5 + 7848iV 4 + 13145iV 3 + 13122iV 2 + 7412iV + 1720 , 
2228iV 10 + 19197iV 9 + 72518iV 8 + 155774iV 7 + 193362iV 6 + 94317iV 5 
-87644iV 4 - 163656iV 3 - 91040iV 2 - 11888iV + 3456 , 
2040iV 15 + 24480iV 14 + 116165iV 13 + 254533iV 12 + 78119iV 11 
-1089300iV 10 - 3414794iV 9 - 5743128iV 8 - 6358562iV 7 - 4824553iV 6 
-2448740iV 5 - 783540iV 4 - 213184iV 3 - 155568iV 2 - 97344iV 
-22464 . 



(30) 



(31) 

(32) 
(33) 



(34) 
(35) 

(36) 
(37) 

(38) 



(39) 



We compared ag'°(JV), a^(N), a^>'\N), a^(N), and a%£(N), Eqs. ( U QQl M M 



-PS,(3),0, 



-PS,(3),0 



NS,(3),0, 



i(3),0, 



[301 . to the fixed moments computed in Ref. [10] and found agreement. 

The OMEs A P ^'q\n) and A^q(N) receive contributions oc rifT F CA,F only. Due to this we 
present as well the constant parts of the renormalized OMEs. They read : 



a 



PS,(3),0 
qq,Q 



(N 2 + N + 2f 



(N — l)iV 2 (l + iV) 2 (2 + N) 
16Rx(N) 



80 



256 



27(iV - l)iV 3 (l + iV) 3 (2 + N) 

R 3 (N) 

~243(iV - l)iV 5 (l + iV) 5 (2 + N) 4 



27 



(S 3 1 +3S 1 S 2 + 2S 3 )+—C 

y 



32R 2 (N) 



(N — l)iV 4 (l + iV) 4 (2 + N) 3 



Si 



(40) 



8 



with 



Ri(N) = 

R2(N) = 

R 3 (N) = 



and 



a 



(3),o 
qg,Q 



with 



AON 7 + 185N 6 + 430iV 5 + 521N 4 + 452iV 3 + 404N 2 - 1QN - 96 , (41) 
233iV 10 + 1744iV 9 + 5937N 8 + 11454iV 7 + 14606iV 6 + 15396N 5 + 1203(W 4 
+3272N 3 - 928iV 2 - 96iV + 288 , 

-42560iV 13 - 445792iV 12 - 2124448N 11 - 6005792iV 10 - 11345024N 9 
-15758592N 8 - 17045248iV 7 - 13567040N 6 - 6545312iV 5 - 1096768iV 4 
+374528iV 3 + 109056iV 2 + 32256iV + 27648 , 



(42) 



(43) 



n f T 2 JC F 



N 2 + N + 2 



Y Si + P 5sfifl + l S2S * + I s * + 2V* + 

LOAT 3 + 13AT 2 + 29 AT + 6) f 3 , 32(5A^ 3 - 16AT 2 + N - 6) 

81iV 2 (l + A0(2 + iV) L^i+^2^iJ+ 81A T2( 1 + A r)( 2 + iV) bs 
^8(109N 4 + 291iV 3 + 478N 2 + 324N + 40) 



N(N + l)(N + 2) 
16(l(W 3 + 13iV 2 + 29iV + 6) 



27iV 2 (l + iV) 2 (2 + Ar) S2 
+ 8(215N 4 + 481iV 3 + 930iV 2 + 748N + 120) 



81iV 2 (l + N) 2 (2 + AT) 
64(A^ 2 + A^ + 2)i? 5 (A^) 



^111 > c 

243A^ 2 (1 + A^) 3 (2 + AT) 1 



9(N - 1)N 3 (1 + A^) 3 (2 + N) 2 2A3(N - 1)A^ 6 (1 + A0 6 (2 + N) 5 

N 2 + N + 2 [ 56 128 160 4 8 2 

O4 — D_4 + —— O3D1 — -D 2 + ~020 1 



+c A 



N(N + l)(N + 2) 



9 



9 



27 



9 



4 o4 64 
— Oi — Jo — 



128 



27 1 9 
32(5N 4 + 20N 3 



S: 



„q 1 

9 ' 9 

t2 



64 256 

— ^2,1,1 - —^3^1 



41iV 2 + 49AT + 20) 3 
liV(l + A0 2 (2 + A0 2 L 1 2,1 " 2 lJ 



64 (5N 4 + 38N 3 + 59N 2 + 31JV + 20) 128 (5N 2 + 8N 
+ 81 N(l + AQ 2 (2 + N) 2 >5 + ~ 1 ' ■ 



10) 



27~Af(l + A0(2 + A0^ 3 
512 (N 2 + N + 1)(N 2 + N + 2) 16R?(N) 
f 9 (AT — 1)A^ 2 (1 + A^) 2 (2 + A^) 2 ~~ 81A/"(1 + AQ 3 (2 + N) 3 2 
_32(121AT 3 + 293AT 2 + 414Ar + 224) _ R 8 (N) 2 

81A/"(1 + A^) 2 (2 + AT) 2 ~ 81AT(1 + AQ 3 (2 + N) 3 1 

16go(A0 „ 8R l0 (N) 

243(A" — 1)A^ 2 (1 + A^) 4 (2 + A^) 4 1 2A3(N - 1)N 5 (1 + N) 5 (2 + NY 



(44) 



R 4 (N) = 24368N 5 + 81984N 4 + 179200N 3 + 225232N 2 + 126880N + 21504 , (45) 
R 5 (N) = 3N 6 + 9N 5 - N 4 - 17N 3 - 38N 2 - 28N - 24 , (46) 
R 6 (N) = 13923N 17 + 180999N 16 + 1064857N 15 + 3812487N 14 + 9348807N 13 

+ 16391845N 12 + 20248499N 11 + 17070917N 10 + 11536274N 9 + 11303496N 8 
+ 13846104N 7 + 16104128N 6 + 22643488N 5 + 29337472N 4 
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+26395008N 3 + 15388416N 2 + 5612544N + 995328 , (47) 
R 7 (N) = 139iV 6 + 1093N 5 + 3438iV 4 + 5776N 3 + 5724iV 2 + 3220iV + 752 , (48) 
R S (N) = 1648N 6 + 11104AT 5 + 34368N 4 + 63856iV 3 + 71904iV 2 + 43264iV + 10880 , (49) 
Rg(N) = + 1244Af 10 + 10557N 9 + 40547N 8 + 90323 + 11449 5N 6 + 49344N 5 

-69902N 4 - 1 15200N 3 - 64352N 2 - 1 1264N + 864 , (50) 
R 10 (N) = 3315N 15 + 39780N 14 + 194011N 13 + 471164N 12 + 416251N 11 - 860568N 10 

-3525799N 9 - 6015120N 8 - 6333994N 7 - 4373672N 6 - 1907512N 5 

-499824N 4 - 217952N 3 - 264192N 2 - 160128N - 34560 . (51) 

In both the constant terms of the renormalized OMEs Eq. (140)440 ( 2 does not contribute anymore. 
Phenomenological applications of the corresponding massive Wilson coefficients are given in 
Ref. [I7j. 

3.2 The Operator Matrix Elements for Transversity 

Transversity is a twist-2 flavor non-singlet operator matrix element related to a tensor oper- 
ator, which cannot be accessed in deep-inelastic scattering, but via polarized semi-inclusive 
deep-inelastic scattering and the polarized Drell-Yan process. The anomalous dimensions for 
transversity are known to NLO [12] and for a series of moments to 3-loop order [31]. Phe- 
nomenological aspects of transversity have been reviewed in Ref. [13] . The moments A" = 1...13 
of the 3-loop massive OME were calculated in Similar to the flavor non-singlet massive 
OME in the vector case we computed the 0(n/) contributions for the transversity operator. The 
constant part of the unrenormalized 3-loop OME is given by 

-tr,(3),o T2r /64 Q 448 32 320 160 



112 640 4 55552 

"C3 + T^T^2 + 77C2 nnn 



9 ' 27 9' 729 
2(3917N 4 + 7834N 3 + 4157N 2 - A8N - 144) 
243N 2 (1 + N) 2 



(52) 



The expression for general values of A" agrees with the corresponding contributions to the mo- 
ments calculated in [IT] before. It is interesting to note that for this color factor the vector and 
tensor operators (I28II52P lead to the same structures in the harmonic sums as for a^'g^ . 



3.3 The Mathematical Structure of the Operator Matrix Elements 

The UfTpCF, A~contributions at 0(a 3 ) to the massive operator matrix elements contain nested 
harmonic sums up to weight w = 4. This also applies to all individual Feynman diagrams, cf. [44J. 
In intermediary results, generalizations of harmonic sums occur, see Appendix |A] As has been 
observed in the computation of various other physical quantities before, such as anomalous 
dimensions and massless Wilson coefficients to 3-loop order [91H9JH5], unpolarized and polar- 
ized massive OMEs to 2-loop order |15j . the polarized and unpolarized Drell-Yan and Higgs- 
boson production cross section, time-like Wilson coefficients, and virtual- and soft corrections 
to Bhabha-scattering [IB] , the classes of contributing harmonic sums are always the same. They 
depend on the loop-order and the topologies of Feynman diagrams involved. 
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In the present case the following harmonic sums emerge 



Si 

S*2, S-2 

S3, <J_3, 5*2,1, 5_2,1 

04, D_4, O31, 5-3,1, 



5-2,2, ^2,1,1, 



5_ 



2,1,1 



(53) 



Note that this class, as for the other processes mentioned above, does not contain the index { — 1}. 
Moreover, we used the algebraic relations between the harmonic sums, cf. [H]. Furthermore, 
structural relations exist between harmonic sums, cf. [23JH7], which reduce the set ( 153]) further. 
Here the sums 



5_ 



2,2, 



S 



3.1 



(54) 



are connected by differential relations w.r.t. their argument iV to other sums of ( 153]) . This is also 
the case for all single harmonic sums S± n , n G N, n > 1, using both the differentiation and 
argument-duplication relation, cf. [29J. Due to this Si represents the class of all single harmonic 
sums. I.e. only the six basic harmonic sums 



Si 

52,1, 5_2,1 

5-3,1) ^2, 1,1) 5_2,11 



(55) 



are needed to represent the 3- loop results for the n/Tj,CV, ^-contributions to the OMEs calculated 
in the present paper. In the final representation we refer to the algebraic basis (153]) and consider 
the basis ( 15 5 p for a later numerical implementation. We sorted the respective expressions keeping 
a rational function in N in front of the harmonic sums ( 153)) and (^-values, like (2 and £3. 

The harmonic sums emerge from the series-expansion of hypergeometric structures like the 
Euler B and T-functions and the Pochhammer-symbols in the (generalized) hypergeometric 
functions pFg^a^e), &i(e); 1) in the dimensional parameter e. This leads to single harmonic 
sums first, which, through summation, turn into (multiple) zeta values [40] and nested harmonic 
sums (29] [30]. The principle steps on the way from single-scale Feynman diagrams to these 
structures have been described in Ref. [23] . 

For phenomenological applications the heavy flavor corrections to the structure functions have 
to be known in x-space. Both the evolution of the parton densities and the Wilson coefficients 
have to be computed at complex values of N. The Mellin-inversion is then performed by a 
numerical contour integral around the singularities of the problem |48j . The analytic continuation 
of the harmonic sums to complex values of N is outlined in Refs. [2"3" ] I4T ] |4"9"] . 



3.4 The OMEs in the Small and Large x Region 

In the small x limit the following leading behaviour of the a^'°, a-gq'^' and afg'Q is obtained 



a 



i(3),0 
l Qg 



PS,(3),0 
Qq 



oc n f T F <C A 



oc —UfTpCF 



18400 448 _ 16 , 

1 C3 H C2 

729 27 S 9 S 



111104 896 h 

+ — c 3 



729 



27 



320 

-w 1 



1 

x 



18540^ 
729 



896 
"27 



<3 



736 
127 



<2 



(56) 
(57) 
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a 



a 



a 



a 



PS,(3),0 

qq,Q 

PS,(3),0 
qq,Q 

NS,(3),0 
qq,Q 

i(3),0 
l gg,Q 



7 (3),0 

l qg,Q 

TR,(3),0 
qq,Q 



oc 



-n f T F C t 



111104 ( 896 



oc rifTpCp 



1024 
16 



729 
47 ' 

~27 +Cs 



320 



l 



1 

x 



oc 



-n,T$C F - In- - 



oc n /T F < C. 



81 \a; 

145408 448 
+ — Cs 



oc n /Tp <^ C 



729 
69472 



27 
512 



64 ' 

Y C2 



68608 896 512 

1 C3 H C2 

729 27 s 27 s 



'OC 



-n f TlC F — In 
27 



729 
1 

;r 



27 



<3 



42688 1024 
+ 



729 



27 



C3 



(58) 
(59) 
(60) 

(61) 
(62) 

(63) 



In case of the singlet and pure-singlet terms the leading behaviour is oc 1/x, while in the 
non-singlet cases it is logarithmic. The small- x asymptotics of d^^' and aj^g turn out 
to be the same. The matrix elements are less singular than the leading terms in the Wilson 
coefficients, cf. [9].[50] . 

In the large x limit one obtains the following leading behaviour, 

oc n f T F (C A -C F )-\n\l 



a 



;(3),0 

l Qg 

NS,(3),0 
qq,Q 



X) 



OC UfTpC F 



27 

55552 448 , 160 ' 

<3 + "7^2 



729 



27 



27 



1 



(1-x). 



~(3),0 (3),0 
a qg,Q' a qg,Q 



4 



oc -n f T 2 F (C A -C F )-\n 4 (l-x) 



TR,(3),0 rp2n 



55552 448 160 
729 ~ ~27 + ~27 



where 



and 



dx- 



a 



a 



a 



PS,(3),0 
Qq 

PS,(3),0 
qq,Q 

PS,(3),0 



[L-X) + JK 

oc —rifTpC F 







1 — X 



oc UfT F C F 
oc njTpC F 



16 5? 
27 iV 2 " 
_ 128 5? _ 



32 



Fi,s(s)-&] 



256 



27 iV 2 
80 5f 



9 

160 



[5i, 3 (x) - C4] 



[5i l3 (^) - C*] 



(64) 
(65) 
(66) 
(67) 

(68) 

(69) 
(70) 
(71) 



•qq,Q — ~r 27 jy2 — Q 

cf. [29j|5l]. In the latter case regular values are obtained for x — Y 1, where 51,3(2) denotes a 
Nielsen integral [52J, 

1 f 1 d? 

S n , p (x) = (-l^^—^ I 7 1" M (^(1-m) • (72) 

The large x limits for fig]' and a^g, resp. ag^ 3 ' 1 ' and d^^'° in the n/ term differ by a factor 

of —8 and —1/8, while the contributions to S^g 3 ^' and o^g 3 '*' are the same. All terms are less 
singular compared to the massless cases [9]. 
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4 The Contributions to the Anomalous Dimensions 



The anomalous dimensions appear in the 1/e term of the unrenormalized OMEs, see Ref. [TO] . 
As all other contributions to this term are known, they can be derived by comparing with the 
1/e terms of the present computation. 



4.1 Vector Operators 

From the OMEs A§J(e, N) and A® Q (e,N) one obtains : 



> 19 



n 2 T 2 



with 



(iV+l)(JV + 2) 



, ^2 ^ ^ I 128 n 32 o3 128 n 64 n 



32 



3iV 



S 2 Si - 



128(5iV 2 + 8iV + 10) 
9N 



64(5iV 4 + 26iV 3 + 47iV 2 + 43iV + 20) 



64(5iV 4 + 20iV 3 + A1N 2 + 49N + 20) , 



9N(N + l)(JV + 2) 
16P 2 (iV) 



9N(N + 1){N + 2) 
64Pi(iV) 



5, 



27iV(iV + l) 2 (iV + 2) 2 



5i 



27(iV - l)N*(N + 1) 3 (N + 2) 3 
32N 2 + N + 2 



+C 



9 



N 



{10S' 3 -5' 1 3 -35' 1 5' 2 } 



32(5iV 2 + 3iV + 2) 32(10iV 3 + 13iV 2 + 29iV + 6) 



-S- 



3N 2 9N 2 
32(47iV 4 + 145iV 3 + 426iV 2 + 412iV + 120) 
27N 2 (N + 1) 



si 



4P 3 (iV) 



27(iV - 1)N 5 (N + l) 4 (iV + 2) 3 



P X (N) = 19iV 6 + 124iV 5 + 492iV 4 + 1153iV 3 + 1362iV 2 + 712iV + 152 , 
P 2 (N) = 165N 12 + 1485N 11 + 5194N 10 + 8534iV 9 + 3557iV 8 - 8899iV 7 

-10364iV 6 + 6800iV 5 + 25896iV 4 + 30864iV 3 + 19904iV 2 

+7296iV + 1152 , 

p 3 (jV) = 99AT 14 + 990iV 13 + 4925iV 12 + 17916A^ n + 46649iV 10 + 72446iV 9 
+32283iV 8 - 95592iV 7 - 267524iV 6 - 479472iV 5 - 586928iV 4 
-455168iV 3 - 269760iV 2 - 122112JV - 27648 . 



(73) 
(74) 

(75) 

(76) 



The n 2 -contribution to the pure-singlet anomalous dimension results from Aq q (e,N) and 

PS,(3' 

gg,Q 



^ p s,(2) 

Iqq 



n}T 2 C F 



32 



(JV- W + l).(iV + 2)l-f (N> + N + mS? + S>) 



+■ 



64 



P 4 (iV) 



-Si 



64 



9 N(l + N)(2 + N) 27 N 2 (l + N) 2 (2 + Nf 



(77) 
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P 4 (jV) = 68iV 5 + 37Ar 6 + 8Ar 7 -lljV 4 -86jV 3 -56jV 2 -104jV-48 , (78) 
P 5 (N) = +52 N 10 + 392 A^ 9 + 1200 A^ 8 + 1353 iV 7 — 317 iV 6 — 1689 A^ 5 

-2103 N 4 - 2672 N 3 - 1496 N 2 - 48 N + 144 . (79) 

Both the 0(rif) contributions to 7^ and 'Jqq'^ have thus been obtained by two independent 
new calculations. 

The n 2 -contribution in the flavor non-singlet case is derived from A^'q\e, N) : 
in f f * ^ g 6 21 27 27 AT 3 (1 + AO 3 J 



with 



p 6 (JV) = 51Af 6 + 153AT 5 + 57Af 4 + 35Af 3 + 96Af 2 + 16Af-24 . 



The anomalous dimensions agree with the moments, resp. the general results, in Refs. [TU1I19| I53]. 
Due to the algebraic compactification we obtain a lower number of harmonic sums Sg(N) if 
compared to Ref. [19], and agree with [35]. For the flavor non-singlet case the anomalous 
dimension has been predicted in [T8] . 

4.2 Tensor Operator 

The contribution to the transversity anomalous dimension oc n 2 is obtained from the single pole 
term of i^ 3) , 

tr,(2) 2 TV f 128 640 128 8 (17lV 2 + 17 N - 8) ] 

The results for the anomalous dimensions constitute a first independent check of the result 
obtained in [TUEl] ■ Again for this color factor the vector- and tensor operators lead to the same 
structures in the harmonic sums. 



5 Conclusions 



We calculated the 0(nf) contributions to the massive operator matrix elements at 0(a 3 ) con- 
tributing to the heavy flavor Wilson coefficients of the deep-inelastic structure function F 2 (x, Q 2 ) 
and to transversity in the asymptotic region for general values of the Mellin variable A^ in the 

PS (3) a (3) 

MS-scheme. Two of the 3-loop OMEs, A qq Q and A qg Q, are known completely now. The Feyn- 
man diagrams contributing are characterized by one massive and (at least) one massless fermion 
line, with both bubble- and ladder-topologies. The local operator insertions are linked to two 
fermion lines and a number of gluon lines. The computation of the Feynman parameter integrals 
has been performed directly by representing the integrals as nested sums over generalized hy- 
pergeometric functions, which result into multiple nested sums over products of hypergeometric 
expressions and harmonic sums. The sums have been solved by applying modern summation 
technologies in difference and product fields. Although in intermediary results in part of the cal- 
culation generalizations of harmonic sums occurred, the final results can be represented in terms 
of rational expressions of the Mellin variable A^ and of harmonic sums of maximal weight w=4. 
The harmonic sums contributing show the same structural pattern as being observed in all other 
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massless 2- and 3-loop calculations. Applying also the structural relations, six harmonic sums 
span the results. The small- and large x behaviour of the constant parts of the OMEs has been 
investigated. In both cases a less singular behaviour than for the massless Wilson coefficients is 

PS (3) (3^ 

observed. The OMEs A q and A y 'q, being completed, do not contain the constant ( 2 after 
renormalization. All results were compared to the fixed moments given in [10] . We mention that 
the present calculation is technically very different from that of computing fixed moments carried 
out previously. From the single pole parts in the dimensional parameter e of the unrenormalized 
OMEs one may derive the respective contributions to the 3-loop anomalous dimensions, which 
are obtained in three cases as a first independent recalculation, using a different method. We 
confirm the results in the literature, both in the deep-inelastic case and for transversity. 
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A Examples for sums occurring in the calculation 



In the present calculation numerous single- to triple finite and infinite sums of an extension of 
the hypergeometric type had to be calculated. For these sums, depending on various summation 
parameters, rij, the ratio of the summands, except the part containing harmonic sums, 

a( 7 n * + 1 - ) , Vi (A.l) 
a(...,n;,...) 

is a rational function in all variables r^. Sums of this type can be represented by basic sums of a 
certain type, which are transcendental to each other and form sum- and product-fields, cf. [25] 
and references therein. The general form of these sums is 



iVi(JV) N 2 (k!,N) N 3 (kuk 2 ,N) 

E E E piSAihfaN)), 

ki=l &2 = 1 &3 = 1 



s al (s 4 (^iv)))r 



s 1 (k i ,N),...,s p (k i ,N) 
s p+1 (ki,N), ...,s p+q (ki,N) 



(A-2) 



with P(x\, x 2 , x%, X4) a polynomial from Q(ki, k 2 , ks, &4, N)[x%, x 2 ,x$, X4], with §i(ki,N), 
. . . , §4(ki, N) and N), . . . , s p+q (ki, N) for some p, q G N being integer linear in ki, k 2 , £3, N, 
with a] an index set, and with the upper bounds Ni(N), N 2 (ki, N), N 3 (ki, k 2 ) being either 00 or 
being integer linear in its arguments. The generalized T-function, cf. [22] , usually includes both 
Beta-functions and binomials. 

In the present calculation one faces more complicated sums than occurring in earlier two-loop 
calculations up to 0(e), [15]. Partly they may reach higher weight than appearing in the final 
result. In the following we present a few examples. 



N-2 00 
ji = l n=l 



N-2\ S 2 (-j l +n + N) 
Ji J n\j x -N-2) 



6 - 23iV + 9iV 2 + 2N 3 
2(N- l) 2 N 2 (l + N)(2 + N) 



+ 



1 



27(-l 



N + 2 (N -1)N(N +1)(N + 2) 



Si 



N(N + 2) 
1 



+ 



48(-l) 



v 



N + 2 (N -1)N(N + 1)(N + 2) 



S3S2 



2S 2 



N(N + 2) 



, ^7(12 + QN -37N 2 + 6N 3 + N 4 ) 
1 1 J 20{-l + N) 2 N 2 {l + N){2 + N) 



1 



,iV 



21 



7 



10(iV + 2) 



5 X + 



10iV(iV + 2) 



< 2 2 + 



+ 



+ 



3(-l) 



TV 



+ 



(N — 1)N(N + 1)(N + 2) N + 2 

3 18(-1) N 

N + 2 (N -1)N(N + 1)(N + 2) 

2 



Si 



S* 



5(-l + N)N(1 + N)(2 + N) 

6 - 23iV + 9iV 2 + 2N 3 
2(-l + iV) 2 iV 2 (l + iV)(2 + N) 

' OA 



N(N + 2) 



2S 2 2 (-l) N (3N-l) 



N + 2 



(N — 1) 3 N 3 



Si 



2 + N 



:SnS 



20-3 



N 3(l2 -6N- UN 2 + 7N 3 + 12N 4 + N 5 ) 
' (-1 + N) 3 N 3 (1 + N)(2 + N) 



16 



+(-1) 



N 



95 2 



+ 



+ 



[-1 + N)N(1 + N)(2 + N) 
3(-l) N 1 



, N/v 3(6-23N + 9N 2 + 2N 3 ) 
+ (- 1 ) 



(N — l)N(N + 1)(N + 2) N + 2 



1)j y / 2 ( 12 - 3 J Ar + 9iv2 + 4Ar3 )^ + { _ 1)N 



-1 + N) 2 N 2 (1 + N)(2 + N) 

5*2 ^5 2 j 

24 



Si 



+ 



25, 



J,2 



(-1 + N) 2 N 2 (1 + N)(2 + N) 



12(-1) 



-1 + + (2 + N) 



Si 



£3,1 



N + 2 



(AT — 1)N(N + l)(iV + 2) iV + 2 



S41 



25_ 2 S_ 2 ,i 45_ 3 ,_ 2 



2 + N AT + 2 



-(-1) 



iV 



-(-!)" 

+(-1)^ 
-(-If 

+(-1)" 

+ "("I 

+(-ir 
(-1) 
+(-ir 



12 



2(6-23A^ + 9AT 2 + 2A^ 3 ) 
-1 + iV) 2 JV 2 (l + N)(2 + N) 

25-2,1,-2 



+ 



N 



+ < 



+ 



If 
("I 



+ 



-1 + N)N(1 + N)(2 + N) 
1 



Si 



52,1,1 - 



N + 2 



-1 + N)N(1 + N)(2 + N) 
30 



425 2 , 2il - 2453,1,1 + 545 2 ,i,i,i 



-1 + N)N(1 + N)(2 + N) 
30 



5 3 5! ( - ) 5! (2) 



-1 + N)N(1 + N)(2 + N) 

N 2(Q-12N + 7N 2 + N 3 ) 
(-1 + JV) 3 JV 3 

2(6-23iV + 9iV 2 + 2iV 3 ) 



5i5i ( - ) 5 3 (2) + 



+ (-1) 



N 



(— 1 + N) 3 N 
6S 2 



(-1 + N)N(1 + N)(2 + N) 



-1 + N) 2 N 2 (1 + N)(2 + N) 
2(6-23iV + 9iV 2 + 2iV 3 ) 



Si 



(-1 + N) 2 N 2 (1 + N)(2 + N) 
6 



5i(2) 
f (-lf- 



12 



-1 + N)N(1 + N)(2 + N) 



Si 



5 2 (2) 



N 



N 



1 + N)N(1 + iV)(2 + N) 

125] 

-1 + N)N(1 + N)(2 + N) 

2(Q - 23N + 9N 2 + 2N 3 ) 
(-1 + N) 2 N 2 (1 + N)(2 + N) 



5 3 (2)^5i,i(-,l 



2(6-23iV + 9iV 2 + 2iV 3 ) ' 
1 ' {-1 + N) 2 N 2 {1 + N){2 + N) 1 



5i 



12 



-1 + iV)iV(l + iV)(2 + A0 
2 (6 — 12 A 7 " + 7N 2 + A 7 " 3 ) 



A? 



(-1 + N) 3 N 3 
2(6-23N + 9N 2 + 2N 3 ) 



-1 + N) 2 N 2 (1 + N)(2 + N) 



Si 
Si 



5i,i(-,l) }S hl (2,l) 



6S 2 



(-1 + A^)A^(1 + AT)(2 + A^) 



17 



+(-1)"; 
-(-!)"- 



30 



-1 + N)N(1 + N)(2 + N) 
6 

-1 + N)N(1 + N)(2 + N)' 



S-2 



Sl,2 



-,1 5i )2 (2,l) 



+ 



x 



, N 2(6-23iV + 9iV 2 + 2iV 3 ) _ ( ^ N 
y ' (-1 + N) 2 N 2 {1 + N)(2 + N) 



S L , U\ 



12 



-1 + N)N(1 + N)(2 + N) 



Si 



■1) 



N 



-1 + N)N(1 + N)(2 + N) 



30S 1 - S li3 (2,l) + 36S M 



2' 



+ 



+ < 



-30S 2 



+ 305,! 



2' 



5 2) i(l,2) 



-1) 



A? 



+ ("!) 



iV 



2(6-23iV + 9A^ 2 + 2iV 3 ) 
(-1 + N) 2 N 2 (1 + N)(2 + N) 

24 



(-1 + A^)A^(1 + A^)(2 + AT) 
("1) N 



Si 



Si 



•&,i(2,l) 



;-l + AT)AT(i + A^)(2 + A^) 



12S M ( -,1 ) S 2 ,i(2,l) + 30S 2 , 3 



-125- 



2,3 



2,^ ) - 185! 



5 3)1 (2,1) + 305 3)2 -,2 -3054,1 



2' 



+ 



2 (6 - 127V + 77V 2 + TV 3 ) _ 
J (-l + iV) 3 AT3 1 iJ 

N 2(6-23iV + 9AT 2 + 2Ar 3 ) 



6S 2 



(-1 + N)N(1 + N)(2 + N) 



-1) 



A? 



(-l + A^) 2 iV 2 (l + A^)(2 + A^) 
30 

J 2 



+ 



i-l) 



(— 1 + N)N(1 + N)(2 + N) 
N 2(Q-23N + 9N 2 + 2N 3 ) 



S 1X1 -,1,2 -,2,1 



;-l + A^) 2 A^ 2 (l + A^)(2 + AT) 



-(-!)"- 
+(-1)^- 



12 



— 1 + N)N(1 + AT) (2 + A?") 
12 

-l + A0Af(l + A0(2 + A0' 



Si 



Si[ j )5i,i,i(1.2,l) 



Si,i -,1 Ui,i(2,l,l) 



+ 



x 



' 2(6-23iV + 9iV 2 + 2iV 3 ) _ 

1 J (-l + AT) 2 Ar 2 (l + A^)(2 + Ar) 1 J i 



12 



-1 + + N)(2 + JV) 



5i 



5i,i >2 [-,1,2 -5 1>1>2 2,-,l -2S 1>li2 2,1,- 



18 



-1) 



N 



[-1 + N)N(1 + N)(2 + N) 



6651,1,3 ( ^,1,2 1 +3651,1,3 I J, 2,1 



+305i,i,3 l,-,2 +305i,i,3 1,2, 



+ 



x 



' , ])N 4(6-23iV + 9iV 2 + 2iV3) _ 
1 1 (-1 + N) 2 N 2 (1 + N)(2 + N) 1 J 



24 



5i, 2 ,i -,2,1 -5 1A1 2,-,l --5 



1 « 



'1,2,1 



2,1, 



(-1 + N)N(1 + N)(2 + N) 
1 



Si 



-(-If 



12 



(-l + N)N(l + N)(2 + N) \2 
1 



5i(-)5i, 2 ,i(l,2,l) 
1 



(-1 + N)N(1 + N)(2 + N) 



30S 1 



2,2 



1,2 +365-12,2 -,2,1 



+4851,3,1 U> 2 > 1 + 30 ^i,3,i 1, n, 2 + 3051,3,1 1, 2, - + 305 2il , 2 -, 2, 1 



+305 2 ,i, 2 1, -, 2 - 125 2 , 1)2 2, -, 1 - 245 2 ,i, 2 2, 1, - + 305 2 , 2 ,i 1,2,- 



245 2 ,2,i 2, -, 1 - 125 2 , 2 ,i 2, 1, - + 3053,1,! -, 1, 2 + 305 3 ,i,i -,2,1 



+ 



(-1)"- 

H-i) N T 



2(6-23iV + 9N 2 + 2N 3 ) 
-1 + N) 2 N 2 (1 + N)(2 + N) 

*1 



■1 + N)N(1 + N)(2 + N) 



1 



5i, 1,1,1 ( -, 1, 2, 1 



251,!,!,! ( 2, -, 1, 1 ) - 251,1,1,1 ( 2, 1, -, 1 ) - 251,1,1,1 ( 2, 1, 1, - 
(-1)" 



(-l + iV)iV(l + A0(2 + A0 



1 



665i,i,i, 2 -,1,2,1 +485!,i,i, 2 -,2,1,1 



+305i,!,i, 2 1, -, 2, 1 + 305!,!,!,, 1, 2, -, 1 + 305!,i, 2 ,i -, 1, 1, 2 



+125^,2,1 ( -,1,2,1 j + 4851,1,2,1 ( -, 2, 1, 1 ) + 305i,i, 2 ,i ( 1, -, 1, 2 
+305i,i, 2 ,i ( 1, 2, 1, - 1 + 3051,2,1,1 (-1,1,2)+ 305i, 2 ,i,i ( |,1,2,1 



+125i, 2 ,i,i -,2,1,1 + 3QS W 1, 1, -, 2 + 305!, 2 ,i,i 1, 1, 2, 



+305 2 ,i,i,i 1, -, 1, 2 + 3052,1,!,! 1, -, 2, 1 + 305 2 ,i,i,i 1, 1, -, 2 



+305 2 ,i,i,i 1,1,2,- + 3052,1,1,1 1, 2, -, 1 + 305 2 ,i,i,i 1, 2, 1, - 



19 



-24S 2>1>1>1 ( 2, 1 1, l) - 2452,1,1,1 (m, ^ l) - 245 2 ,i,i,i (2, 1, 1, i 
-125i,i,i,i,i ( i 1, 2, 1, l] - 365i,i,i,i,i (\, 2, 1, 1, 1 



+ (-1) 



A' 



35? 



(-l + iV)iV(l + JV)(2 + iV) 



+ (-1) 



N 



2 + 9N - 5N 2 



[-1 + N) 2 N(1 + N)(2 + N) 



+ (-1) 



A' 



+ 



6-11N + 2N 2 
-1 + N) 2 N 2 (2 + N) 

3(-l) N 



Si 



+ 



1 



(N -1)N(N + 1)(N + 2) N + 2 
95 2 



S2K3 



+C2 (-1) 



N 



2(-l + N)N(1 + N)(2 + N) 



+ (-1) 



N 2 + 3N- 2 2+N N -2N 2 -3- 2 1+N N 2 + 6N 3 - 2 1+N N 3 - 3N 4 



(-1 + N) 3 N 2 {1 + N){2 + N) 
-12 + N + 27N 2 - AN 3 



2(-l + N) 2 N 2 (1 + N)(2 + N) 



+ 



6(-l) 



A? 



AT + 2 (7V-l)7V(7V + l)(7V + 2) 



Si 



1 



iV(iV + 2) 
25_ 2 



'(So 



6(-l) 



A? 



iV(iV + 2) 
25_ 2 



N + 2 
'ivT2 



5i + 



-1) 

35_3 
N + 2 



(N -1)N(N + 1)(N + 2) iV + 2 



N -6 + 3N + 18N 2 - 20N 3 - 3iV 4 + 2N 5 
(-1 + N) 3 N 3 (1 + N)(2 + N) 

12 



+ (-1) 



N 



+ 



"("I) 



A? 



(-1 + N)N(1 + N)(2 + N) 
3S 2 



5 2 ,i 



■1 + N)N(1 + N)(2 + N) 



, , N 6- 12N + 7N 2 + N 3 . . N -Q + 23N -9N 2 -2N 3 
+(-1) r~, An o An — + (-1) 



+ 



(-1) 
(-1)" 



+ 



(-l + iV) 3 iV 3 v ' (-1 + N) 2 N 2 (1 + N)(2 + N) 

N -6 + 23iV - 9N 2 - 2iV 3 
(-1 + N) 2 N 2 (1 + N)(2 + N) 

6 c~ 

{-\ + N)N{\ + N){2 + N) 1 

N 6 - 23N + 9N 2 + 2N 3 

(-1 + iV) 2 iV 2 (l + N){2 + N) 



Si 



5i(2) 



5 2 (2) 



+(-1) 



N 



6 



-5i 



(-1 + N)N(1 + N)(2 + N) 

(-1 + JV)iV(l + AT) (2 + N) 
(-If 



5i,i(2,l) 



-1 + N)N(1 + N)(2 + N) 



35i, 2 (2, 1) - 155 2 ,i(l, 2) + 65 2 ,i(2, 1) - 65 1 , 1 , 1 (2, 1, 1) 
20 



with 



+ 



18(-1) 



A? 



+ 



N + 2 (N — 1)N(N + i)(N + 2) 
27(-l) N 9 



(N — 1)N(N + 1)(JV + 2) 2(iV + 2) 



5(a,6) 



(iV - l)iV(iV + l)(N + 2) Sl [ 2 ' 5l ' 1,2(2 ' X ' 1} ' 



T(a + 6) 



(A.3) 



(A.4) 



(-i)^so»(-^" 2 )^i(i)^iW 



AT-2 -j+N-2 oo 

^ S^+" n )^ +n+1 )^ +n + 2 )^ +n + 3 )^ 1 " iV - 2 ) 



(-1) 



N -1 + 2N-3N 2 , -16 + 12A^ + 10A^ 2 -17A^ 3 -31A^ 4 



+ 



(-1)" 

N 



N 3 (l + N)(2 + N) 
1 



+ 



8N 3 (1 + N)(2 + N) 



1 



N(2 + N) (N + 1)(JV + 2)J 4(N + 2) 



+ 



1 



s 4 



2(N + 2) 



(1 + N)(2 + N) 



S-9 + So 



-8-2N + N 2 + 5N 3 + 5N 4 
4AT 3 (1 + A^)(2 + A^) 



+ 



+ 



(-1) N (N~1) 
N 2 (N + 2) ' 2N 2 (N + 1) 

2(-l) N 



Si + (-i) 



N 



N(2 + N) 



S. 



+ (-1) 



N 



1 + 2N 2 \ 2 + N + N 2 

(N + l)(N + 2) ( ' v ^ N 2 (l + N)(2 + N) ~ 3 + 1 [2N 2 (1 + N)(2 + N) 



+ (-1) 



N 



2(-l + N) 
N 2 (2 + N) 



N{2 + N) 



S. 



3(~l) iV c 
+ iV(2 + A0 4 



+ 
+ 



-1) N (1 -N) 



+ 



N 2 (N + 2) 2(N + 2) 

S 3 ,i + 



S2,l - ("I) 



N 



N(2 + N) 



S 



2,-2 



N(N + 2) 2(N + 2) 
-(-1) 



(-1) 



„ 2(-l + 2iV 2 ) 



N 2 (l + N)(2 + N) 



5-2,1 + (-1) 



iV(2 + A") 

at r 2 



Si 



N(2 + AO 



-3,1 



+ 



N(2 + N) S ~ 2 ' hl 
N 3 (N + 1)(A" + 2)J 



16 + 4N - 2N 2 + A^ 3 + A" 4 
8A^ 3 (1 + AT)(2 + AT) 



Si 



2A^ 2 (A^+ 1) 



C 2 , 



(A.5) 



:-l)*jBV,n)(-t + ?-*)S 2 (n) 



N-2 -j+N~2 oo 

5 S 5 (j+n)(j+n + l)(j+n + 2)(j+n + 3)(j 1 -iV-2) 



5| 



A(N + 2) 



(-l) N S 2 
2N(N + 2) 



+ 



:-i) 



iV 



AT(AT + 2) 2AT2(AT+1) 



Si 



21 



iV 2 



'2N 2 {N+1) 



+ (-1) 



N 



N(2 + N) 



+ (-1) 



N 



8 - 19iV + 24N 2 



(-1 + N)N 2 (1 + N)(2 + N) 



S-0 + 



+ 



1 



So 



2N 2 (N+1) 
+(-l) N 
+Si 



-1) 



N 



2N(N + 2) 4(^ + 2) 
1 



(N + l)(N + 2). 

-48 - 24iV + 71N 2 + 95iV 3 
48AT 2 (l + iV)(2 + A^) 

S, 



(1 + N)(2 + N) 
2 



N(2 + N) 

I (-1)" 
N(N + 2) 

H-i) N 



N(2 + N) 



S. 



S. 



S-2 + 



1 



N(2 + N) 



N 



N(2 + N) 



■S. 



SfS-2 + 



N(2 + N) 



S. 



S 2 ,i + 

1 



-(-1) 

' 2{-l) N 
N(N + 2) 



:Si + 



N(2 + N) 2N 2 (N+1) 



-1) 



A? 



N(2 + N) 



Si 



S-2,1 



S-,i + 



+ 



+ 



+ 



(-1) 



N 



N(2 + N)~ 3,1 ' N(2 + N) 
—2 + N — 2N 2 



52,1,1 + 



iV(2 + N) 



-2,1,1 



2 + N - N 2 - 2N 3 

+ 7TT + 



Si 



2(-l + N)N 2 {1 + N){2 + N) 2N 2 {1 + N)(2 + N) 2N 2 (N + 1) 



(-1) 



N 



So 



-1) 



JV 



+ 



iV(2 + N) 



S. 



C-2 



2(N + 2) N(N + 2) 

' -12- QN + N 2 + N 3 
12N 2 (1 + N)(2 + N) ^ (N — 1)N 2 (N + 1)(N + 2) 



(A.6) 



j=l j 1 = l 



A? 



(j + 2)(ji — N — 2) 
1 



+ 



2(JV + l)(JV + 2) 2(iV + 2) 
-3 - 3iV - N 2 



S 2 



Sf 



2(N + 2) 



+ (-1) 



jv — 4 — 5A" — 3iV 2 — A^ 3 



+ (-1) 



(l + iV) 2 (2 + TV) 2 ' V v iV(l + N) 2 (2 + A^) 2 

+ 28N + 37N 2 - A2N 4 - 38N 5 - UN 6 - 2N 7 



Si 



N ' 



2N 2 {1 + N) 3 {2 + N) 3 
+ 28N + A9N 2 + 39N 3 + 6N 4 - ION 5 - 6N 6 - N 7 



+ (-1) 



N 



N 2 {l + N) 3 {2 + N) 3 
2 (-8-28N -27N 2 -8N 3 )S 2 



(1 + A0(2 + A0 



S-2S2 + 



2Ar 2 (l + A^)3(2 + AO 3 



4 + 5N + 3N 2 + N 3 + 1 ^6 + 2A^-8A^ 2 -6A^ 3 - A^ 4 



A^(l + A r ) 3 (2 + AT) 1 



(1 + A0 3 (2 + AT) 3 
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+(-1) 



N 



-4-3N + 2N + 3N 3 + N 4 
N(l + N) 2 (2 + N) 2 

3 



+ 



(1 + N)(2 + N) 



Si 



S* 



2(-A + 2N 2 + N 3 ) 



2(1 + N)(2 + N) * N 2 (2 + N) 2 

3( - 4 - 3N + 2N 2 + 3N 3 + iV 4 ) 



N(l + N) 2 (2 + N)' 



S, 



-1 



N 2(A + 5N + 3N 2 + N 3 ) 



N(l + N) 2 (2 + Nf 



S- 2 + (-1 



16 + SN - AN 2 - N 3 
N 3 (2 + N) 3 

S-3 



1 



(N + l) 2 (N + 2) 



+ (-1 



,iv 



+ 



(1 + JV)(2 + JV) 
2 4 

5-4 + 



[1 + N)(2 + N) 



[1 + N)(2 + N) 



S 



2,-2 



+ 



2(-A-3N + 2N 2 + 3N 3 + N 4 ) 
N(l + N) 2 (2 + N) 2 

6 4 



;i + ao(2 + ao 



5-3,1 + 



1 + JV)(2 + iV) 



(1 + JV)(2 + JV) 

5-2,1,1 



Si 



s. 



2.1 



(A.7) 



In the above examples also so-called generalized harmonic sums occur [3~T|I32]. They obey 
the following recursive definition : 



S mit ...(xi, • ••; N) 



^a^a,.-^) ^; 

i 1= l 1 i 2 =l 2 

+5 mi -i- ?Tt2i?Tt3i ... (xi • x 2 ,x 3 , AT) . 



(A.8) 



The sums 5* may be reduced to nested harmonic sums for Xi G { — 1,1}. In the present calcula- 
tion the values of extend to {—1/2, 1/2, —2,2}. These sums occur in ladder like structures, 
cf. 0,[35], but may also emerge if contributions to 3-loop Feynman diagrams, containing a 2- 
point insertion, are separated into various terms. They were even observed in case of the more 
complicated massive 2-loop graphs [15] if large expressions are arbitrarily separated. In part of 
the sums terms ex 2 N , which lead to an exponential growth in the large N limit, occur. However, 
all these contributions cancel for each individual diagram. In the present case the weight of these 
sums can reach w = 5 intermediary, depending on the ^-structure of the contribution, although 
only w = 4 sums will emerge in the final results. Examples for these sums are : 



Si(l/2,iV), S 2 (-2;N), S 2> i(-1,2;N), S 3 ,i( 
Si,i,i )2 (-1, 1/2, 2, -1; iV), S 2j3 (-2, -1/2; N), 
S 2l 2,i(-1,-1/2,2;A0, etc. 



■1/2; AT), 



(A.9) 



The algebraic and structural relations for these sums are worked out in Ref. [32J. Similar to 
the case of harmonic sums, corresponding basis representations are obtained. They allow to 
simplify involved structures as of Eq. (jA.30 and finally lead to the reduction of the results for 
the individual diagrams to a representation just in terms of nested harmonic sums. The nested 
sums emerging in this work, which were not given before in Refs. [15J and those being closer 
related to the structure of harmonic sums [30], are of the type illustrated above. The latter 
have been calculated using C. Schneider's packages Sigma [25], Evaluatelnf initeSums [54J and 
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J. Ablinger's package HarmonicSums |2"5] . 
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